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ABSTRACT 

We analyse the environment of low redshift, z < 0.2, SDSS quasars using the 
spectral and photometric information of galaxies from the Sloan Digital Sky Survey 
Third Data Release (SDSS-DR3). We compare quasar neighbourhoods with field and 
high density environments through an analysis on samples of typical galaxies and 
groups. 

We compute the surrounding surface number density of galaxies finding that 
quasar environments systematically avoid high density regions. Their mean environ- 
ments correspond to galaxy density enhancements similar to those of typical galaxies. 

We have also explored several galaxy properties in these environments, such as 
spectral types, specific star formation rates, concentration indexes, colours and active 
nuclei activity. We find a higher relative fraction of blue galaxies in quasar envi- 
ronments compared to groups and typical galaxy neighbourhoods. Consistent with 
this picture, the distribution of the concentration index of these galaxies also indi- 
cate a larger fraction of late-type objects. By analysing the available information of 
galaxy spectra we have also studied the distribution of the star formation rates of 
these neighbour galaxies finding that quasar environments are populated by objects 
with an enhanced star formation activity. An analysis of the relative flux ratios of 
[OIII]\5700/H(3 and [NII]X6583/Ha of emission-line galaxies shows no an excess of 
nuclei activity in quasar neighborhood with respect to the environment of a typical 
galaxy. 

We conclude that low redshift quasar neighbourhoods (r p < l/i _1 Mpc, AV < 
500kms _1 ) are populated by bluer and more intense star forming galaxies of disk- 
type morphology than galaxies in groups and in the field. Although star formation 
activity is thought to be significantly triggered by interactions, we find that quasar 
fueling may not require the presence of a close companion galaxy (r p < 100/i _1 kpc, 
AV < SSOkms" 1 ). 

As a test of the unified AGN model, we have performed a similar analysis to 
the neighbours of a sample of active galaxies. The results indicate that these neigh- 
bourhoods are comparable to those of quasars giving further support to this unified 
scenario. 

Key words: quasars: active galaxies : statistics- distribution - galaxies: groups: 
general - 



their neighbourhoods in the local Universe. It is well known 
that several properties of galaxies depend on the environ- 
ment where they formed and evolved where a variety of 
processes such as star formation, tidal stripping, merging, 
etc. can determine the nature of galaxies. It should also be 
considered the possibility of AGN feedback processes which 



1 INTRODUCTION 

Important clues on galaxy formation and evolution may 
be obtained by characterizing statistically the properties of 
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could induce significant changes in the evolution of their 
companion galaxies (see for instance Croton et al. 2005). 

Previous works aimed to characterize the quasar neigh- 
bourhood indicate that high redshift quasars can be used as 
signspots to search for rich density regions (Djorgovski 1999, 
Hall & Green 1998, Fukugita 2003). However, at lower red- 
shifts (z < 0.3) different studies indicate that quasars could 
reside in environments similar to those of normal galaxies 
(Smith, Boyle & Maddox 1995, Sorrentino et al. 2006) or 
they could be located in groups (Fisher et al. 1996) or clus- 
ters of galaxies (Mclure & Dunlop 2001). Coldwell et al. 
2002, analysing both the projected cross-correlation func- 
tion and colours of galaxies around a sample of quasars 
and AGNs found that their typical galaxy density environ- 
ment corresponds to groups of galaxies. Moreover, Coldwell 
& Lambas (2003, hereafter Paperl), using the spectral in- 
formation of the 2dF Galaxy Redshift Survey (2dFGRS), 
found that the galaxies within r p < l/i _1 Mpc from quasars 
have a stronger star formation activity. These results are 
also supported by findings of Soechting et al. (2002, 2004) 
who found that low redshift quasars follow the large-scale 
structure traced by galaxy clusters but they are not placed 
in the central area of galaxy clusters. More likely, they are 
in the cluster periphery or between two, possibly merging, 
galaxy clusters. 

The large recent spectroscopic surveys, the Sloan Dig- 
ital Sky Survey and the 2dF Galaxy Redshift Survey allow 
us investigate a wide range of galaxy properties such as mor- 
phology, spectral types, colours and, also, the galaxy den- 
sity using volume-limited samples. In this paper we extend 
the work of Paperl by analysing the environment of Sloan 
quasars using the spectroscopic and photometric informa- 
tion available for the Sloan Digital Sky Survey Third Data 
Release (SDSS). We divided the sample in two range of red- 
shift 0.03 < z < 0.1 and 0.1 < z < 0.2 in order to anaelyse 
different effects in these two redshift ranges populated by 
galaxies of different luminosities. 

The layout of the paper is as follows. In section 2 we 
describe the data, section 3 shows the analysis of the local 
density estimates, section 4 and 5 describes the statistical 
analysis performed with the photometric and spectroscopic 
data respectively, and in section 6 we provide a brief discus- 
sion of the main results. 



Mj = -23, PSF magnitude i < 19.1 and highly reliable 
redshifts (Schneider et al. 2003). 

The galaxy groups from SDSS were identified by 
Merchan & Zandivarez (2005) by using the friends-of-friends 
algorithm developed by Huchra & Geller (1982) which was 
improved by implementing a procedure to avoid the arti- 
ficial merging of small systems in high density regions and 
applying an iterative method to recompute the group centres 
position. The group sample has a median velocity dispersion 
of 230kms~ and we restrict our target to those groups with 
a minimum number of 8 members to obtain higher density 
environment. 

As a way to reject the hypothesis that quasars could 
reside in environments such as that corresponding to galaxy 
groups we use three different target samples, taken from 
SDSS, in order to compare quasar environment with those 
corresponding to typical galaxies and galaxy groups and we 
make appropriate comparisons of galaxy characteristics of 
quasar neighbourhoods with respect to those in a low density 
environments corresponding to typical galaxies and denser 
environments as galaxy groups. 

We analyzed galaxy properties in the neighbourhood of 
different target samples in two ranges of redshifts in order 
to explore for a luminosity dependence: 

0.02 < 2 < 0.1, hereafter Z\: 418 SDSS Quasars, 1147 
SDSS Galaxies and 779 SDSS Galaxy Groups . 

0.1 < z < 0.2, hereafter Z2: 1652 SDSS Quasars, 1153 
SDSS Galaxies, 102 Galaxy Groups. 

The targets were selected to match the observed quasar 
redshift distribution, showed in Fig. 1, in order to have unbi- 
ased and directly comparable results and the samples have 
the largest number of objects (within the redshift restric- 
tion) which suffices to provide reliable statistical results. By 
doing so, we assume that there are not redshift-dependent 
systematics such as those associated to low luminosity galax- 
ies. A Kolmogorov-Smirnov test yields that the redshift dis- 
tributions are very similar with a 90 % level of confidence. 

Tracer galaxies consist of all objects within projected 
distance r p < 3ft" 1 Mpc and with radial velocity difference 
AV < 500kms~ 1 relative to the targets, so that both targets 
and tracers have a similar redshift distribution. 



2 DATA 

The Sloan Digital Sky Survey (SDSS) in 5 optical bands 
will map one-quarter of the entire sky and perform a redshift 
survey of galaxies, quasars and stars. The third data release, 
DR3, from SDSS provides a database of 374767 galaxies and 
51027 quasars with measured spectra. The five filters u, g, 
r, i, and z cover the entire wavelength range of the CCD 
response (Fukugita et al. 1996). The main galaxy sample is 
essentially a magnitude limited spectroscopic sample (with 
a Petrosian magnitude) ru m < 17.77, most of galaxies span 
a redshift range < z < 0.25 with a median redshift of 0.1 
(Strauss et al. 2002). 

The quasar sample is defined by quasars which have 
at least one emission line with a full width at half max- 
imum larger than lOOOkms -1 , luminosities brighter than 



3 LOCAL DENSITY ESTIMATES 

A useful characterization of the local galaxy density can 
be obtained by measuring the distance to the N th near- 
est neighbour and estimating the density within that dis- 
tance. The advantage of this method is to use a systemat- 
ically larger scale in lower-density regions which improves 
sensitivity and precision at low densities. This is a two- 
dimensional estimate but we use the redshift information 
to lower the projection effects. We choose a fixed velocity 
interval of AV = lOOOkms -1 to compute the local density 
which correspond to galaxies within ~ 3a from the centre 
of a galaxy cluster (Balogh et al. 2004) and this allows in- 
clusion of galaxies in system with large velocity dispersion. 
When computing projected distances we have assumed a 
flat cosmology (go = 0.5, Q = 0.3, fix = 0.7) and a Hubble 
constant Hq = 100/ikms _1 Mpc _1 . 
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to typical galaxies. The Ei density parameter corresponds to 
a closest neighbor distance estimate (within the luminosity 
restrictions) so that the results of Fig. 2 indicate that galaxy 
interactions are not likely to be directly associated to quasar 
phenomena. The errorbar in this figure and in all figures 
were calculated by using boostrap error resampling (Barrow 
et al. 1984). The results shown in Fig. 2 strongly suggest 
that local quasars avoid systematically high and moderately 
high density regions such as groups of galaxies. 




Figure 1. Rcdshift distributions of targets in the redshift inter- 
vals Zl and Z2. Dashed lines correspond to quasars, solid lines 
to typical galaxies, and dotted lines to groups of galaxies. 




Figure 2. Distributions of logiQ Si (left panels) and logio 
£5 (right panels) for galaxies in the neighbourhoods of quasars 
(dashed lines), typical galaxies (solid), and groups of galaxies 
(dotted), for two ranges of rcdshift. 



We calculate surface densities, Ei and E5, correspond- 
ing to the projected distance d\ and ds of the first and fifth 
neighbour brighter than M r < —20.5 respectively. This ab- 
solute magnitude limit assures completeness within the red- 
shift range explored, z < 0.2 
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In Fig. 2 we show the relative distribution of Ei and 
E5 where it can be appreciated that the local density of 
quasar environment is pretty similar to that corresponding 



4 PHOTOMETRIC PROPERTIES 

Different magnitude measurements are provided in SDSS, 
asinh magnitude, Petrosian magnitude, fiber magnitude, etc. 
where galaxies bright enough to be included in the spectro- 
scopic sample, r < 17.7, have relatively high signal-to-noise 
ratio. Since Petrosian magnitudes are model independent 
and yield a large fraction of the total flux, roughly con- 
stant with redshift, they provide an adequate magnitude 
measurement. We have used the modified form of the Pet- 
rosian (1976) system which measures galaxy fluxes within a 
circular aperture whose radius is defined by the shape of the 
azimuthally averaged light profile in order to measure a con- 
stant fraction of the total light, independent of the position 
and distance of the objects. 

4.1 SDSS Colours 

Galaxy colours can be used as estimators of the galaxy evo- 
lution. In clusters, the large fraction of red galaxies indicate 
an old population of galaxies with a low star formation rate. 
Galaxies in poor groups or in the field are bluer and with 
stronger star formation rate. 

We calculate the corrected colours it — r of galaxies in 
the neighbourhood of the targets by using K and extinction 
corrected absolute magnitude (Blanton 2003). We analyse 
the relative distribution of u — r for galaxies around the 
targets where we observe an excess of blue galaxies in the 
environments of quasars. To determine the significance of 
this excess we calculate the relative fraction of galaxies bluer 
than a given threshold value, we selected the fraction of 
galaxies with u — r < 2.3 which correspond approximately 
to the mean u — r for the spectroscopic survey at z < 0.2. 
The results are shown in Fig. 3. 

We can see in this figure that the fraction of blue neigh- 
bours of a typical galaxy tends to decrease at small sep- 
arations. This effect that can be interpreted as the usual 
morphology-density relation and, as expected, this decrease 
is even stronger as one approach group centres. However, 
contrary to this behaviour, the blue fraction does not de- 
cline at close separations and there is a significant relative 
excess of blue galaxies in quasar neighbourhoods compared 
to typical galaxies. This result is perhaps an indication of 
feedback process operated by quasars which could delay star 
formation at early stages allowing for more recent activity 
from the remaining gas. 

4.2 Morphology 

In SDSS are also available the apparent radii containing 50% 
and 90% of the Petrosian flux for each band. The ratio of 
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Figure 3. Left panels: Fraction of blue galaxies, u — r < 2.3, 
relative to the total number in each bin of projected distance r p 
from the targets: quasar neighbourhoods (empty circles), typi- 
cal galaxies (solid lines), and groups of galaxies (filled triangles). 
Right panels: Distributions of u — r colours of neighbour galaxies, 
within r p < 0.5ft -1 Mpc, to Sloan quasars (dashed lines), groups 
of galaxies (dotted lines) and typical galaxies (solid line) in the 
range 0.02 < z < 0.1, b) 0.1 < z < 0.2. 



Figure 4. Left panels: Fraction of disk-type galaxies, C < 2.5, 
relative to the total number in each bin of projected distance r p 
from the targets: quasar neighbourhoods (empty circles), typi- 
cal galaxies (solid lines), and groups of galaxies (filled triangles). 
Right panels: Distributions of concentration index C of neighbour 
galaxies, within r p < 0.5h -1 Mpc, to Sloan quasars (dashed lines), 
groups of galaxies (dotted lines) and typical galaxies (solid line) 
in the ranges Zl and Z2. 



these fluxes, the concentration index C, is correlated with 
morphology. Since galaxies with a de-Vaucouleurs profile 
have a value of C ~ 3.3 and disk galaxies have a concentra- 
tion index C ~ 2.4 this parameter can be used as a simple 
morphological classifier. In Fig. 4 we can see the distribu- 
tion of the concentration index parameter and the fraction 
of late type galaxies defined with C < 2.5 as a function of r p . 
The larger relative fraction of disk-type galaxies in quasar 
neighbourhoods can be clearly appreciated corresponding to 
~ 20% excess of disk-type objetcs. 



5 SPECTRAL PROPERTIES 
5.1 Spectral Type Classification 

The regulation of star formation in galaxies can be strongly 
influenced by close companions. Besides the effects of en- 
vironment, galaxies in pairs have enhanced star formation 
over a control sample with similar characteristics which is 
stronger for galaxy pairs in field that in groups of galaxies. 
Lambas et al. (2003) and Alonso et al (2004) found clear ev- 
idence that star formation is enhanced over 40% in close in- 
teractions. This star formation induction increases for small 
relative velocity and projected separation r p . 

On the other hand, feeding a massive black hole at the 
centre of galaxies may require particular conditions which 
could be strongly influenced by environment. Although at 
large redshifts quasars and radio galaxies are frequently as- 
sociated to clusters, at lower redshifts, a dense environment 
may be hostile to the presence of massive black holes (Cold- 
well, Martinez & Lambas 2002, Coldwell & Lambas 2003). 
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Figure 5. Fraction of emission-line galaxies (s > 0.2), Fi, and 
fraction of early galaxies (s < 0), F2, relative to the total num- 
ber in each bin of projected distance r p from the targets in the 
neighbourhoods of quasars (empty circles), typical galaxies (solid 
line) and groups of galaxies (filled triangles) in the ranges Zl (left 
panels) and Z2 (right panels). 



Properties of galaxies in SDSS Quasar environments at z < 0.2 5 



In Paperl we explored the nature of galaxies in the 
vicinity of the different target samples using the 2dFGRS 
spectral type index r/. This parameter approximately delin- 
eates the transition between early and late morphological 
types for r\ ~ —1.4 . When considering r\ > 3.5 we are deal- 
ing with galaxies particularly active, such as starbursts with 
recent episodes of star-formation or AGNs. 

SDSS galaxies are also classified by a principal compo- 
nent analysis, PC A (Connolly & Szalay 1999), where five 
eigencoefficients, ecoeff, are extracted. Similarly to the 2dF- 
GRS spectral type classification SDSS provides a spectral 
type parameter s with the first two eigencoefficients like 
s = arctan(—ecoef f2/ecoef fl) which ranges from about 
-0.35 to 0.5 for early to late galaxies. Taking in to account 
this and the corresponding galaxies for SDSS and 2dFGRS, 
we consider galaxies with s < as being similar spectra 
than their 2dF counterparts with (77 < —1.4). Also, SDSS 
galaxies with s > 0.2 correspond to the 2dF galaxies with 
77 > 1.1 The fraction, Fi, galaxies with strong star formation 
activity (s > 0.2), and the fraction F2 of early spectral type 
galaxies (s < 0), are shown for the two ranges of redshifts Zl 
and Z2 in Fig. 5. The effect is similar to that of 2dF galax- 
ies in Fig. 1 of Paperl where it is clear that galaxies around 
quasars differ from that in galaxy and group environments. 

5.2 Star formation vs. AGN activity 

Several important properties of galaxies have been derived 
for subsamples of SDSS: stellar masses, indicators of recent 
major starbursts, current total and specific star-formation 
rates (Brinchmann et al. 2004) and emission-line fluxes 
(Tremonti et al. 2004) both for the regions with spectroscopy 
and for the galaxies as a whole; gas-phase metallicities; AGN 
classifications (Kauffmann et al. 2003) based on the stan- 
dard emission line ratio diagnostic diagrams, etc. 

The relation between spectral lines, [Oil I] A5007, HP, 
[A//JA6583 and Ha luminosities, can be used to analyse 
possible dependence of the relative numbers of AGNs and 
star forming galaxies in quasar neighborhoods compared to 
galaxy and group environments. We constructed the Bald- 
win, Phillips & Terlevich (BPT, 1981) line-ratio diagram and 
we used the Kauffmann et al. (2003) criteria to differentiate 
AGNs galaxies from other emission-line objects dominated 
by star formation. An AGN is defined if 

log([OIII]/H/3) > 0.61/(log([NII/Ha\) - 0.05) + 1.3 (2) 

Figure 6 shows the BPT diagram for redshift range Zl 
where it can be appreciated a lack of differences in the distri- 
bution of AGN neighbours for quasars, galaxies and AGNs. 
Similar results were found for the redshift range Z2. 

To quantify this result we calculated the percentage 
of AGN neighbours within l/i -1 Mpc from the three target 
samples obtaining that 37.9 % , 35.2% and 34.5 % of galaxies 
in the vicinity of groups, galaxies and quasars, respectively, 
are AGNs indicating that the AGN activity is not strongly 
affected by quasars. 

Taking into account the previous results which indicate 
that quasar environments are populated by blue, disk-type, 
star forming galaxies, we have computed the fraction of 
AGNs (in different environments) which have these char- 
acteristics (g — r < 0.7 and C < 2.5). The corresponding 
fractions of AGNs obtained are 6.7 % , 5.3% and 6.4 % in 




Log(NII/Ha) 

Figure 6. The distribution of the neighbour (r p < l/i _1 Mpc) 
galaxies in the BPT line-ratio diagram for the three target sam- 
ples at Zl. The solid line shows the division of the samples. 

the vicinity of groups, galaxies and quasars, respectively in- 
dicating that AGN contribution to blue star forming galaxies 
is not significant in any environment. 

5.3 Star Formation Rates 

The methods for deriving the star formation rates are based 
on models where the contribution of the nebular emission by 
HII regions and diffuse ionized gas combined and described 
in terms of effective metallicity, ionization parameter, dust 
attenuation at 5500 A, and dust to metal ratio (Bruzual & 
Chariot 1993, Chariot et al. 2002). Taking into account these 
issues, Brinchmann et al. (2004) provide accurate total star 
formation rates estimates free from aperture bias. Moreover, 
Kauffmann et al. (2003) developed a method to constrain 
star formation histories, dust attenuation and stellar masses 
of galaxies based on two stellar absorption lines indices, the 
4000 A break strength and the Balmer absorption line index 
Hcca. 

Our main interest here is to analyse the logarithmic 
specific star formation rate logSFR/M* [log yr _1 ], where 
M* is the estimated mass in star, for galaxies in the differ- 
ent target environments. In the right panels of Fig. 7 it is 
shown the distribution of logSFR/M* from which it can be 
seen that quasar environments have a population of galaxies 
with a higher star formation rate than field galaxies (galaxy 
groups show a low rate of star formation as expected from 
the systematic decline of SFR with local density). 

In order to quantify the excess of star formation activity, 
we calculate the fraction of galaxies with log(SFR/M*) > 
— 10.0 and the results are shown in Fig. 7 (left panels). Sev- 
eral theories have proposed that galaxy-galaxy interactions 
fuel the AGN activity by driving gas into the cores of galax- 
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Figure 7. Left Panels: Fraction of galaxies with strong specific 
star formation rate, logSFR/M* > —10.0, relative to the total 
number in each bin of projected distance r p from the targets in 
the neighbourhoods of quasars (empty circles), typical galaxies 
(solid line) and groups of galaxies (filled triangles) Right Panels: 
Distribution of star formation rate per mass unit, logSFR/M* , 
within r p < 0.5h _1 Mpc from the targets in the neighbourhoods 
of quasars (dashed line), typical galaxies (solid line) and groups 
of galaxies (dotted line) . 



ies and thus onto the black holes. We have directly explored 
whether quasars have companions more frequently than the 
control samples. Taking into account the observed thresh- 
olds in projected separation r p and relative velocity AV 
for galaxy interactions to effectively induce star formation 
we calculated the fraction of targets with close companions 
(r p < lOO/i^kpc and AV < SSOkms" 1 ). We find a low 
fraction (< 15%) of close companions associated to quasars, 
similar to that observed for galaxies in general. This result 
indicates that quasar fueling may not require the presence of 
close companion galaxies (consistent with the distribution of 
Ei shown in Fig. 2). Moreover we have calculated the distri- 
butions of star formation excluding these close neighbours 
(those within r p < lOO/i^kpc and AV" < 350kms _1 ). The 
results are very similar to those of Fig. 7, indicating the lack 
of relevance of close interactions in driving these relations. 



5.4 Comparison of Quasar and AGN 
environments. 

The unification hypothesis for active galactic nuclei and 
quasars (Antonucci 1993) is a model where all active galaxy 
classes are fundamentally the same phenomenon seen at 
different orientations which are not due to intrinsic differ- 
ences As a final analysis we have explored the characteristics 
(colours, specific star formation rates and concentration in- 
dexes) of galaxies in AGN environments (ie. taking AGNs as 
targets) in a similar fashion as in previous sections serving 
as a test of the unified AGN model. The AGN target sam- 
ples were taken from Kauffmann et al. (2003) restricted to 





r„ (Mpc h-') (Z,) 



r p (Mpc h- 1 ) (Z a ) 



Figure 8. Comparison of quasar and AGN environments. Left 
panels correspond to the rcdshift interval Zl and right panels to 
Z2. a), b) and c) correspond to the fractions of blue colour index 
galaxies (u — r < 2.3), disk-types (C < 2.5), and strong specific 
star formation rate objects (logSFR/M* > —10) respectively. 
The empty circles correspond to quasar environments and the 
stars to AGN environments. The solid lines show the results of 
typical galaxy neighbourhoods. 



have a similar redshift distribution than the quasar, group, 
and typical galaxy samples analysed previously. The results 
are shown in Fig. 8 where it can be appreciated the simi- 
larity of galaxy properties surrounding quasars and AGNs. 
The fact that the environments of quasars and AGNs are 
comparable provide further support to the unified model. 
We notice, however, that AGN environments more closely 
resemble that of quasars in the high redshift range, a fact 
that may be related to the higher luminosity of AGNs here. 



6 SUMMARY 

We have performed a statistical analysis of local quasar en- 
vironments using SDSS data, our main conclusions are: 

• Nearby quasars systematically avoid high density 
regions. Local surface density estimates show that SDSS 
quasars reside in similar density enhancements to typical 
galaxies, significantly less dense than cluster environments. 

• Star formation activity in the surrounding of quasars 
is higher than in the neighbourhood of typical galaxies. This 
important property may provide a link between star forma- 
tion and the onset of quasar activity. 

• Quasar environment is populated by galaxies system- 
atically bluer, and with a disk-type morphology. 

• In spite of the fact that galaxy interactions efficiently 
trigger star formation, we find no statistical evidence that 
the presence of close companions is associated to quasars. 



Properties of galaxies in SDSS Quasar environments at z < 0.2 7 



•The presence of quasars does not affect the fraction of 
AGNs in neighbour galaxies. 

• The results mentioned are present in the two redshift 
ranges explored z < 0.1 and 0.1 < z < 0.2 although they 
are stronger at lower redshifts, indicating that the brightest 
galaxies are not likely to be the major contributors to the 
effects. Also, there might be a signal for evolution in the 
sense that the effects are stronger locally. 

• AGN and quasar environment are similar regarding 
the galaxy characteristics explored and agree best in the 
highest redshift interval giving support to the unified sce- 
nario of AGNs and quasars. 

Although quasar hosts tend to be red and bulge- 
dominated galaxies (Kauffmann et al. 2003), the indication 
of a relative excess of quasar neighbour galaxies with a large 
gas fraction could suggest either an external origin for the 
accretion onto the central black hole, or that quasars may 
affect substantially the evolution of star formation up to the 
Mpc scale. 
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